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ABSTRACT 


Faraday  rotation  data  from  Satellite  S-66  were  accumulated  at  the 
Electronics  Laboratory  of  the  Military  Research  and  Development  Center, 
Bangkok,  Thailand  from  November  1964  through  June  1967,  a  period  of 
intermediate  sunspot  number  (sunspot  minimum  condition  at  1964/1965  and 
sunspot  maximum  condition  at  1968/1969).  These  data  provide  a  means  of 
improving  ionospheric  frequency  predictions  for  regions  in  Thailand 
remote  from  Bangkok  (i.e.,  where  no  ionosonde  data  are  available). 

In  this  report,  two  methods  of  analysis — involving  rotation  rate 
and  total  number  of  rotations — are  applied  to  the  Bangkok  observations. 

The  rotation-rate  method  is  used  to  determine  the  electron  content  when 
the  angle  between  the  ray  path  and  geomagnetic  field  is  90°  (transverse 
position).  The  total-rotations  method  is  used  to  determine  latitudinal 
variations  of  electron  content  of  the  equatorial  ionosphere  for  the 
selected  passes  from  a  joint  analysis  of  observations  at  three  stations 
in  Thailand — Songkhla,  Bangkok,  and  Chiengmai. 

The  diurnal  variation  of  the  local  electron  content  shows  low  early- 
morning  values  and  high  afternoon  values,  resulting  in  a  large  maximum- 
to-minimum  ratio.  A  secondary  maximum  was  observed  during  the  nighttime. 
The  diurnal  and  seasonal  variations  of  electron  content  are  examined 
together  with  the  variation  with  solar  activity  as  measured  by  the  10.7-cm 
solar  radio  flux. 

The  diurnal  development  and  collapse  of  the  equatorial  anomaly  in 
electron  content  is  clearly  seen  from  the  total-rotations  data.  The 
Chapman  equation  is  used  for  conversion  of  these  total  electron  content 
values  to  peak  electron  density  in  order  to  compute  the  F-layer  critical 
frequency  (foF2)  and  its  variation  as  a  continuous  function  of  latitude. 
The  computed  foF2  values  are  compared  with  those  obtained  from  three 
ionosondes  which  were  located  at  Songkhla,  Bangkok,  and  Chiengmai,  and 
the  agreement  is  typically  within  10  percent. 
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FOREWORD 


The  work  described  in  this  report  was  performed  with  the  support, 
and  using  the  facilities,  of  the  Military  Research  and  Development  Center 
(MRDC)  in  Bangkok,  Thailand.  The  MRDC  is  a  joint  Thai-U.S.  organization 
established  to  conduct  research  and  development  work  in  the  tropical 
environment.  The  overall  direction  of  the  U.S.  portion  of  the  MRDC  has 
been  assigned  to  the  Advanced  Research  Projects  Agency  (ARPA)  of  the 
U.S.  Department  of  Defense  who,  in  1962,  asked  the  U.S.  Army  Electronics 
Command  (USAECOM)  and  the  Stanford  Research  Institute  (SRI)  to  establish 
an  electronics  laboratory  in  Thailand  to  facilitate  the  study  of  radio 
communications  in  the  tropics  and  related  topics.  The  MRDC-Electronics 
Laboratory  (MRDC-EL)  began  operation  in  1963  [under  Contract  DA  36-039 
AMC-00040(E) ]  and  since  that  time  the  ARPA  has  actively  monitored  and 
directed  the  efforts  of  USAECOM  and  SRI.  In  Bangkok,  this  function  is 
carried  out  by  the  ARPA  Research  and  Development  Field  Unit  (RDFU-T) . 

The  cooperation  of  the  Thai  Ministry  of  Defense  and  the  Thailand  and 
CONUS  representatives  of  the  ARPA  and  USAECOM  made  possible  the  work 
presented  in  this  report. 
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I  INTRODUCTION 


One  of  t>e  major  activities  of  communication  research  in  equatorial 
areas  is  the  study  of  the  ionosphere.  A  recent  phase  of  this  effort  is 
the  study  of  the  development  and  collapse  of  the  equate  jal  anomalies  in 
electron  content  and  F-layer  critical  frequency — the  equatorial  anomaly 
is  the  name  given  to  a  region  of  unusually  high  value  of  electron  concen¬ 
tration  or  critical  frequency  in  the  F-region  of  the  ionosphere  near  the 
magnetic  equator.  To  increase  understanding  of  these  ionospheric  char¬ 
acteristics,  records  of  the  Faraday  rotation  on  40-  and  41-MHz  signals 
from  the  S-66  (Explorer  22)  radio  beacon  satellite  have  been  accumulated 
from  November  1964  through  June  1967  at  the  MRDC  Electronics  Laboratory1 
in  Bangkok,  and  at  two  remote  sites  located  north  and  south  of  Bangkok 
during  1966  and  1967.  The  projection  of  a  typical  satellite  path  of 
descending  pass  on  the  earth  over  these  sites  is  shown  in  Fig.  1.  Data 
were  collected  for  sites  in  the  range  of  geographic  latitude  between 
about  7° N  and  19° N. 

The  methods  of  analyzing  the  data  on  Faraday  rotation  of  radio  signals 
from  a  satellite  are  described  in  detail  in  a  Special  Technical  Report, 

STR  14, 6  which  reported  the  results  of  electron  content  over  Bangkok  and 
the  equatorial  anomaly  during  the  period  of  sunspot  minimum. 

This  report  presents  the  results  of  electron-content  measurements  over 
Bangkok  for  a  period  of  intermediate  sunspot  number  (sunspot  minimum  condi¬ 
tion  at  1964/1965  and  maximum  condition  at  1968/1969).  The  content  is  used 
to  calculate  the  equivalent  slab  thickness  and  the  mean  scale  height.  In 
this  report  the  diurnal  development  and  collapse  of  the  equatorial  anoma¬ 
lies  in  electron  content  and  F-layer  critical  frequency  are  analyzed  from 
selected  Faraday  rotation  data  from  satellite  passes  observed  by  the  re¬ 
ceiving  stations  Songkhla,  Bangkok,  and  Chiengmai.  Using  the  Chapman 
equation,  the  electron  co.itent  data  are  used  to  calculate  the  F-layer 
critical  frequency  of  the  equatorial  ionosphere.  The  assumptions  involved 
are  given,  and  the  possible  sources  of  error  due  to  those  approximations 
are  discussed. 

+ 

References  are  listed  at  the  end  of  the  report. 


II  THEORETICAL  BACKGROUND 


When  a  linearly  polarized  electromagnetic  wave  is  propagated  through 
a  low-loss,  ionized  medium,  such  as  the  ionosphere,  and  is  under  the 
influence  of  the  earth's  magnetic  field,  the  polarization  of  the  resultant 
electric  field  vector  rotates  gradually  as  the  wave  moves  through  the 
ionosphere.  This  phenomenon  is  commonly  referred  to  as  Faraday  rotation. 

The  relationships  between  a  satellite  and  a  ground  observation 
point  are  shown  in  a  simplified  form  in  Fig.  2.  The  important  angular 


FIG.  2  FARADAY  ROTATION  GEOMETRY 


parameters  are  Z,  the  satellite  zenith  angle  (the  angle  between  the  signal 
ray  path  and  the  vertical),  and  9,  the  angle  between  the  signal  ray  path 
and  the  earth's  magnetic  field.  Two  positions  of  the  satellite  are  indi¬ 
cated —  the  usual  case,  where  9  ^  9 (f  ,  and  a  special  case,  where  9  =  9(]P  . 
The  region  near  9  ca  &  is  called  the  region  of  quasi-longitudinal  (QL) 


3 


propagation,  and  the  point  where  6  =  0°  is  called  longitudinal  case.  The 
region  near  0  s»  90°  is  called  the  region  of  quasi-transverse  (QT)  propa¬ 
gation,  and  the  point  where  0  =  90°  is  called  the  transverse  case  or 
position.  (At  this  point  the  sense  of  rotation  will  be  reversed.) 

A  linearly  polarized  signal  can  be  treated  as  the  sum  of  two  circu¬ 
larly  polarized  waves  with  opposite  directions  of  rotation.  It  can  be 
shown  that  the  angle  of  the  plane  of  polarization  rotates  as  the  wave 
travels  through  the  ionosphere. 


\  =  Radio  wavelength  in  the  medium 

ke  =  nekQ  =  Ordinary  (left-hand  polarization) 
component  phase  constant 

kr  =  nrkQ  =  Extraordinary  (right-hand  polarization) 
component  phase  constant. 


From  the  geometry  illustrated  in  Fig.  3,  propagation  is  in  the  z 
direction  end  an  element  of  path  length  in  this  direction  is  d l.  The 
change  in  angle  of  polarization  rotation,  dH,  as  the  wave  proceeds  a 
distance  Al ,  is  given  by: 

k 

dfl  =  g"  (ne  “  nr>  (D 

where  the  resultant  electric  field  of  the  wave  is 

E  =  E  e^^ez  +  E  e^^rz 
®  r 
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In  the  quasi-longitudinal  case,  the  Appleton-Hartree  equation  for 
the  refractive  Index  under  the  QL  approximation  is 


n  1 

e,r 


and  , 

k 

dH  =  (XYL)d£ 


(2) 

(3) 


Z 


In  terms  of  mks  units, 

dQ  =  £-  NH  cos  9  di  (4) 

f 2  ° 

where 

Hq  =  Earth's  magnetic  field  intensity,  amp- turn/m 

3 

N  =  Electron  density,  electron/m 

f  =  Frequency  of  satellite  transmission,  Hertz 

-2 

K  =  A  constant  =  2.97  X  10  ,  mks  units 

0  =  Angle  of  polarization  rotation,  radian 

d l  =  Element  of  path  length 
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0  =  Angle  between  the  magnetic  field  and  wave  normal  (direction  of 

propagation) . 

To  find  the  electron  content  in  a  vertical  column  (vertical  columnar 
electron  density)  through  the  ionosphere,  it  is  convenient  to  replace  the 
path  element  d£  by  dh  sec  x>  where  x  is  the  local  zenith  angle  and  dh  is 
an  element  height.  With  this  substitution,  Eq.  (4)  becomes 

dO  =  NH  cos  9  sec  x  dh  .  (5) 

r  ° 

Integrating  over  height  from  zero  to  the  satellite  height, 
h 

K  »  S 

0  =  —  NH  cos  9  sec  x  dh  .  (6) 

.2  <J  o  o  * 


The  total  number  of  rotations  of  the  plane  of  polarization  of  the 
linearly  polarized  resultant  is 

n  _ 

h 


R  = 


2tt 


K 


H  cos  9  sec  x 


I 


N  dh 


2nf 

Q 


W  . -  n  S 

=  —  B  cos  9  sec  x  ]  Ndh 


(7) 

(8) 


where 


-13 

Q  =  A  constant  =  3.77  X  10  ,  mks  units 

f  =  Frequency,  MHz 


B  =  Earth's  magnetic  intensity,  Gauss 

B  cos0secx  =  A  weighted  mean  over  the  range  of  integration; 
or,  more  compactly, 

R  =  G  J  Ndh  .  (9) 

In  order  to  determine  the  electron  content  using  Eq.  (9),  the  value 
* 

of  G  may  be  calculated  at  any  point  in  the  spatial  coverage  around  the 

observer.  In  general,  the  G  value  must  be  the  weighted  average  of  the 

actual  value  of  B  cos  9  sec  x  along  the  line  of  the  ray  path.  By  using 

the  average  ionospheric  profile,  the  value  of  G  is  determined  near  the 
* 

The  calculation  of  G  values  is  given  in  Appendix  A. 
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height  of  the  centroid  of  the  ionospheric  electron-density  profile.  The 
centroid  height  for  the  area  of  interest  is  obtained  by  fitting  the 
Faraday  rotation  electron  content  data  to  the  ionospheric  distribution 
profile.  Therefore,  the  value  of  G  is  evaluated  at  a  level  approximately 

a 

50  km  above  the  height  of  maximum  density,  which  can  be  determined  from 
a  true-height  analysis  of  vertical-incidence  ionograms  obtained  near  the 
location  and  time  of  the  satellite  passages. 

Equation  (9)  was  derived  under  the  following  assumptions: 

(1)  The  wave  frequency  is  so  high  that  refraction  can  be 
neglected.  (That  is,  the  wave  frequency  is  much  higher 
than  the  maximum  plasma  and  collision  frequencies.) 

(2)  The  effect  of  horizontal  electron  gradients  can  be 
neglected . 

(3)  The  quasi-longitudinal  approximation  holds.7 

(4)  The  vertical  velocity  component  of  the  satellite  can 
be  neglected. 

(5)  The  satellite  antenna  has  a  fixed  spatial  orientation. 

(6)  A  single  ray  path  exists  for  both  modes. 

Equation  (9)  is  a  good  approximation  to  a  more  exact  formula  as 
long  as  the  quasi-longitudinal  approximation  is  valid.  The  quasi¬ 
longitudinal  approximation  is  not  valid  near  9  =  9CP  ,  the  transverse 
position  (which  occurs  at  the  time  Tq ) .  At  9  =  9CP  ,  both  R  and  G  in 
Eq.  (9)  go  to  zero  and  the  expression  for  electron  content  becomes 
indeterminate.  It  is  realized  that  at  the  transverse  point  the 
rotation  is  a  complicated  series  of  polarization  changes.  The  method 
used  to  obtain  the  electron  content  at  the  transverse  position  is  a 
method  of  interpolation — the  rotation-rate  method.  To  obtain  the 
rotation  rate,  Eq.  (9)  is  differentiated  with  respect  to  time,  giving 

g  =  rhs  Ndh .  £  +  o .  s-  r”s  Ndh  .  no 

dt  «  dt  dt  J 

o  o 

When  the  horizontal  gradients  in  the  ionosphere  are  neglected  for 
a  very  narrow  latitudinal  range  of  the  ionosphere,  the  second  term 
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becomes  zero.  Therefore,  the  rotation  rate  becomes  proportional  to  the 
electron  content 


2£  »  22.  '"s  Ndh 

dt  dt  O 


n  • 

P  S  Ndh  =  £  .  (12) 

J  °  o 

In  Eq.  (12)  the  dot  implies  a  time  derivative.  The  G  can  be  determined 
indirectly  by  differenting  the  G  with  respect  to  time,  to  obtain 


d9s  ‘  dt 


where  9  is  satellite  latitude  and  — - 
s  dt 


is  the  time  derivative  of  satellite 


latitude.  Near  the  geographic  equator,  for  Satellite  S-66  at  a  height  of 
d9 

1000  kilometers.  s  has  the  value  0.0565  degrees  per  second, 
dt 

Equations  (9)  and  (12)  are  the  basic  relationships  for  calculating 
electron  content  from  experimentally  obtained  Faraday  rotation  data. 
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Ill  METHOD  OF  ANALYSIS 


A .  GENERAL  APPROACH 

One  widely  used  method  of  studying  the  electron  distribution  in  the 
ionosphere  is  the  Faraday  rotation  technique.  The  principal  advantage 
of  the  Faraday  rotation  technique  for  calculating  electron  content  used 
in  this  study  is  the  simplicity  of  data  acquisition  and  analysis.  Analy¬ 
sis  of  Faraday  rotation  data  can  be  done  in  two  ways — by  the  rotation- 
rate  method  and  by  the  total-rotations  method.  The  rotation-rate  method 
is  suitable  for  analyzing  a  large  number  of  records  and  has  the  advantage 
of  being  relatively  insensitive  to  horizontal  gradients  of  electron 
content.  However,  it  can  be  used  only  in  calculations  of  electron 
content  near  the  transverse  position  and  so  can  give  values  at  only  one 
location.  The  rotation-rate  method  has  been  applied  to  Bangkok  data. 

The  total-rotations  method  is  more  complex  but  has  the  advantage  of 
yielding  values  of  electron  content  over  a  range  of  latitudes.  This 
method  is  easy  to  use  for  data  from  a  station  near  the  magnetic  equator, 
because  the  usual  polarization-rotation  ambiguity  can  be  resolved  when 

the  satellite  passes  through  the  T  position.  By  counting  the  number 

o 

of  rotations  of  the  plane  of  polarization  from  this  transverse  point, 
one  is  able  to  determine  an  unambiguous  measurement  of  electron  content 
for  latitudinal  variation.  The  total-rotations  method  was  applied  to 
the  data  from  the  three-station  chain — Songkhl a-Bangkok-Chiengmai . 

Figure  4  gives  an  example  of  the  appearance  of  the  Faraday  polari¬ 
zation  fading-rate  anomaly  on  signal-amplitude  records  which  were 
observed  at  three  stations  during  a  typical  passage  of  Satellite  S-66. 
These  Faraday  fading  data  were  recorded  on  40  MHz  for  the  purpose  of 
determining  the  latitudinal  variation  of  electron  content  over  Thailand. 
The  fading  pattern  seen  in  Fig.  4  is  caused  by  the  resultant  polarization 
vector  of  the  received  signal  alternately  falling  parallel  to  and  orthog¬ 
onal  to  the  linearly  polarized  antenna.  Since  the  signal  polarization 
(electric  field  intensity)  vector  becomes  parallel  with  the  dipole  twice 
during  each  rotation  of  36CP ,  a  complete  rotation  corresponds  to  two 
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periods  of  the  fading  pattern.  The  unique  point — the  time  of  occurrence 

of  transverse  propagation  — is  plainly  visible  in  Fig.  4.  The  change 

in  the  pattern  near  T  is  caused  by  the  reversal  in  the  direction  of  the 

0 

polarization  vector  rotation.  Local  time  was  recorded  on  the  chart  to 
relate  the  observed  data  to  the  physical  position  of  the  satellite  in  its 
orbit.  The  satellite  position  and  height  at  time  corresponding  to  the 
local  time  were  obtained  from  ephemeris  data  provided  by  the  U.S.  National 
Aeronautics  and  Space  Administration  (NASA). 

B.  ROTATION-RATE  METHOD 

The  rotation-rate  method  of  analysis  is  based  on  the  rate  of  change 

of  rotation  near  the  T  position  to  determine  the  electron  content  at 

o 

that  position.  In  this  method  of  analysis,  Eq.  (12)  is  used  to  obtain 
the  electron  content.  In  this  equation,  R  can  be  measured  directly  from 
the  data  by  graphical  calculation.  The  G  coefficient  is  determined  by 
differentiating  G  as  shown  by  Eq.  (14).  The  calculation  has  been  pre¬ 
pared  for  computer  calculations. 

A  typical  plot  of  the  number  of  rotations  observed  on  a  record  is 

shown  in  Fig.  5.  The  number  of  rotations  (left-hand  scale)  is  counted, 

starting  from  an  arbitrary  first  null  of  the  signal  amplitude  in  the 

Faraday  rotation  record,  and  the  time  corresponding  to  the  third,  fifth, 

etc. ,  nulls  are  used  to  plot  the  number  of  rotations  as  a  function  of 

local  time.  Figure  5  illustrates  the  calculation  of  the  rotation  rate, 

R,  by  a  three-point  numerical  differentiation  method  (using  the  slope  of 

the  curve  near  T  ).  The  number  of  rotations  is  scaled  from  a  1-minute 
o 

period  centered  at  the  T  position.  This  gives  the  rotation  rate  (R) 

o 

directly  in  rotations  per  minute.  The  ratio  of  the  rotation  rate  (R) 
to  the  coefficient  (G)  yields  the  electron  content  as  given  in  Eq.  (12). 


* 

Local  time  (GMT  +  7  hours)  was  established  using  a  General  Radio  115B 
Frequency  Standard  adjusted  using  phase  referenced  to  the  VLF  trans¬ 
mission  from  GBR  in  England  on  15  kHz.  Time  ticks  generated  by  a 
GR  1123A  Syncronometer  from  this  standard  were  compared  daily  with 
time  transmission  from  WWVH  and  JJY  on  10  and  15  MHz.  The  resulting 
time  estimate  is  accurate  to  5  milliseconds. 
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FIG.  5  SAMPLE  ANALYSIS  CALCULATION 


C.  TOTAL-ROTATIONS  METHOD 

For  a  station  near  the  magnetic  equator,  ambiguity  in  the  total 

number  of  rotations  can  be  resolved  easily  for  most  records.  There  is 

a  time  during  each  satellite  passage  when  transverse  propagation  occurs 

and  the  resulting  reversal  of  wave  polarization  produces  a  visible 

perturbation  of  the  Faraday  fading  pattern  at  that  time  (T  ).  The 

o 

total  number  of  rotations  is  defined  as  equivalent  to  zero  at  the  T 

o 

position.  The  total-rotations  method  of  analysis  uses  the  total  number 

of  rotations  from  the  transverse  position  to  a  given  satellite  position 

to  determine  the  electron  content  at  that  position.  The  right-hand 

scale  of  Fig.  5  illustrates  the  total  number  of  rotations  referred  to 

T  for  different  times  during  a  satellite  pass.  The  total  number  of 
o 

rotations  is  counted  from  T  to  the  time  corresponding  to  the  satellite 

o 

latitude  on  both  sides  of  T  .  Satellite  S-66  has  a  nearly  polar  orbit; 

o 

hence  the  satellite  latitude  is  a  more  sensitive  variable  in  this 
analysis . 

To  follow  the  above  procedure,  the  local  time  corresponding  to  the 
satellite  position  must  be  known  accurately.  The  satellite  ephemeris 
is  used  to  tabulate  the  different  satellite  positions  corresponding  to 
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local  time  during  the  observation  period.  The  total  number  of  rotations 
(R)  is  counted  from  the  right-hand  scale  in  Fig.  5  and  the  corresponding 
satellite  position  is  obtained  from  the  local-time  tabulations  of  the 
satellite  ephemeris.  In  this  method  of  analysis,  Eq.  (9)  is  used  to 
obtain  the  latitudinal  electron  content.  A  computer  program  has  been 
prepared  for  this  analysis  that  can  be  used  to  compute  the  latitudinal 
variation  of  electron  content  for  any  satellite  passage  of  interest 
(see  Appendix  B) . 


IV  RESULTS 


A.  LOCAL  ELECTRON  CONTENT 

The  local  electron  content  over  Bangkok  was  calculated  by  using 
Eq.  (12)  at  the  transverse  position.  This  position  corresponds  to  the 
ionospheric  point  of  approximately  14.2°N  latitude  for  all  satellite 
passes  recorded.  The  electron  content  values  presented  in  this  section 
have  been  calculated  from  the  Faraday  rotation  data  obtained  from 
Satellite  S-66  from  November  1964  through  June  1967,  a  period  of 
intermediate  sunspot  number  (sunspot  minimum  condition  at  1964/1965 
and  maximum  condition  at  1968/1969). 

To  see  the  main  features,  the  electron  content  was  plotted  for  the 
ascending  and  descending  passes  separately,  as  shown  in  Figs.  6  and  7 
respectively.  Because  the  satellite  precessed  westward,  the  time  of  the 
satellite  passage  over  Bangkok  became  progressively  earlier.  No  differ¬ 
ence  is  apparent  between  the  results  for  the  two  directions  of  passage. 
The  time  of  day  advances  from  right  to  left  in  Figs.  6  and  7,  and  each 
calculated  value  of  electron  content  is  associated  with  a  particular 
hour  of  the  day.  Therefore,  the  passage  time  moved  through  a  24-hour 
period  in  about  five  and  a  half  months. 

In  addition  to  the  calculated  values  of  electron  content,  F-layer 

peak  electron  density  and  the  daily  mean  10.7-cm  solar  flux  values  are 

shown  in  Figs.  6  and  7.  In  these  figures,  the  top  curve  is  the  electron 

content,  the  middle  curve  is  the  F-layer  peak  electron  density  calculated 

from  the  ionosonde  observations  over  Bangkok,  and  bottom  curve  shows  the 

daily  mean  10.7-cm  solar  flux.  The  main  feature  of  the  electron  content 

results  is  the  large  diurnal  variation — low  early  morning  and  high 

afternoon  values.  The  maximum  value  of  electron  content  varied  from 
16  2 

20  X  10  electron/m  at  the  low  sunspot  number  (February  1965)  to 
16  2 

55  X  10  electron/m  at  high  sunspot  number  (March-April  1967)  and 

occurred  around  1400  hours  local  time,  while  the  minimum  value  of 

16  2 

electron  content  varied  from  1  X  10  electrons/m  at  the  low  sunspot 
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DIURNAL  VARIATION  OF  ELECTRON  CONTENT  AND  F'LAYER  PEAK  ELECTRON  DENSITY  AND  SOLAR 
FLUX  —  DESCENDING  PASSES 


16  2 

number  (January-February  1965),  to  10  X  10  electron/m  at  the  high 
sunspot  number  (March  1967),  and  occurred  in  the  early  morning  hours. 

The  slope  of  the  diurnal  curve  rose  sharply  in  the  morning  hours, 
reached  a  peak  in  the  afternoon,  and  decayed  rapidly  to  a  minimum  in 
the  early  morning  hours.  Secondary  maxima  were  observed  at  the  night¬ 
time  hours.  This  nighttime  behavior  of  ionospheric  electron  content 
agrees  with  the  observations  by  other  investigators  at  other  lati¬ 
tudes.9-12 

The  values  of  the  F-layer  peak  electron  density  (see  Figs.  6  and  7) 
were  calculated  from  ionosonde  data  obtained  at  times  very  nearly  corre¬ 
sponding  to  the  satellite  passes.  The  curves  of  electron  content  and 
peak  electron  density  are  generally  somewhat  similar.  It  is  noticeable 
that  the  secondary  maximum  of  peak  electron  density  also  appeared  at 
nighttime.  The  cause  of  these  secondary  maximum  values  is  not  clear. 

Solar  flux  values  (10.7  cm)  are  plotted  in  Figs.  6  and  7  as  an 
index  of  solar  ultraviolet  activity  during  the  period  of  satellite  (and 
vertical-incidence-sounder)  ionospheric  observations.  The  main  features 
are  the  diurnal  variations  of  solar  flux,  which  show  a  correlation  with 
the  envelope  of  the  peak  electron  content  values, 

B.  SEASONAL  VARIATIONS 

It  is  apparent  in  Figs.  6  and  7  that  the  diurnal  variations  shown 
in  these  plots  are  contaminated  by  seasonal  variations.  The  seasonal 
variations  of  the  daytime  maximum  values  of  electron  content  are  illus¬ 
trated  in  Fig.  8.  In  order  to  study  this  variation  further,  the  daytime 
maximum  electron  content  values  during  sunspot  number  minimum  in  1964 
are  taken  from  Ref.  6. 

The  principal  features  are  the  low  electron  content  value  during 
the  sunspot  minimum  in  1964,  and  the  electron  content  increasing  as 
the  sunspot  number  increases  toward  maximum.  The  seasonal  variation  in 
Fig.  8  shows  the  maximum  values  of  electron  content  occurring  near  the 
March-April  and  September-October  equinoctial  periods,  since  in  this 
period  the  sun  is  directly  overhead  at  the  equator,  and  minimum  values 
occurring  near  June-July  and  December- January  of  solstitial  periods. 
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FIG.  8  SEASONAL  VARIATION  OF  PEAK  ELECTRON  CONTENT 


C.  SOLAR  VARIATIONS 

To  investigate  the  dependence  of  electron  content  upon  solar  activity, 
the  values  of  electron  content  during  the  middle  of  the  day  (Figs.  6  and 
7),  when  the  diurnal  rate  of  change  is  small,  are  plotted  against  the 
corresponding  daily  mean  10.7-cm  solar  flux  values  (see  Fig.  9).  In 
Fig.  9  the  lower  points  are  for  the  low  sunspot  number  period  and  the 
upper  points  are  for  the  high  sunspot  number  period.  This  indicates  an 
increase  of  electron  content  with  solar  flux  that  is  the  same  order  of 
magnitude  as  the  effect  of  sunspot  number.  Apparently  the  electron 
content  is  linearly  dependent  upon  solar  flux,  at  least  for  small-to- 
moderate  values  of  solar  flux. 

D.  SLAB  THICKNESS  AND  SCALE  HEIGHT 

The  slab  thickness  and  mean  scale  height  of  the  ionosphere  are 
closely  related  to  the  electron  distribution  profile.  They  can  be 
computed  from  the  following  equations:1 3 >14 

Ndh  =  4.13  H  (N  F)  (15) 

J  max 
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FIG.  9  SOLAR  VARIATION  OF  PEAK  ELECTRON  CONTENT 

where  h 

f  SH(h)  dh 

»  =  -^h - 

s 

H(h)  =  Scale  height 

H  =  Mean  scale  height 

h  =  Satellite  height 

s 

N  F  =  Peak  electron  content  of  the  F-layer. 
max 

By  definition,  the  total  electron  content  is  equal  to  the  product 
of  equivalent  slab  thickness  and  peak  F-layer  electron  content.  There¬ 
fore,  the  equivalent  slab  thickness  can  be  identified  ss  4.13  H  in  Eq. 
(15). 

It  is  evident  that  the  mean  scale  height  and  slab  thickness  are 

related  by  a  factor  of  4.13.  The  individual  value  of  electron  content, 

determined  by  the  rotation-rate  analysis,  can  be  divided  by  the  peak 

electron  density  (N  F)  obtained  from  ionosonde  data  taken  at  the  same 

max 

time  to  obtain  the  equivalent  slab  thickness.  The  peak  electron  density 
can  be  obtained  from  ionospheric  data  collected  at  Bangkok  by  a  modified 
C-2  ionosonde,  by  using  the  following  relation: 
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(16) 


2  10 

N  F  =  1 .24  (f oF2 )  X  10 
max 

where  foF2  is  the  observed  F  -layer  ordinary  wave  critical  frequency 

‘  3 

(megahertz),  and  N  F  is  the  peak  electron  density  (elcctron/m  ). 

'  max 

The  diurnal  and  seasonal  variation  of  slab  thickness  and  mean  scale 
height  are  clearly  shown  in  Fig.  10.  Data  for  the  ascending  and  de¬ 
scending  passes  are  plotted  separately.  In  Fig.  10  a  definite  variation 
between  day  and  night  is  discernible,  since  both  slab  thickness  and  mean 
scale  height  have  the  diurnal  variation.  The  maximum  daytime  slab  thick¬ 
ness  and  mean  scale-height  values  vary  between  about  300  to  400  km  and 
75  to  100  km  respectively.  The  minimum  nighttime  values  of  slab  thick¬ 
ness  and  mean  scale  height  range  approximately  from  120  to  190  km  and 
30  to  45  km,  respectively.  These  values  are  generally  comparable  to 
other  published  values. 

E.  ACCURACY  OF  RESULTS 

The  rotation-rate  method  is  suitable  for  the  analysis  of  a  large 

number  of  recordings  and  is  relatively  insensitive  to  horizontal  electron 

gradients  near  the  transverse  (T  )  position.  It  is  obvious  from  Eq.  (9) 

o 

that  a  Faraday  effect  cannot  be  claimed  to  exist  at  the  transverse 
position.  The  only  method  to  obtain  the  electron  content  at  the  trans¬ 
verse  position  is  the  rotation-rate  method.  In  the  complete  equation, 

Eq.  (10),  the  second  term  (which  is  neglected  for  the  rotation-rate 
analysis)  increases  the  rotation  rate  when  horizontal  gradients  are 
present.  By  using  the  horizontal  gradients  from  the  three-stations 
analysis  of  the  selected  passes,  we  have  determined  that  the  second 
term  causes  less  than  5-percent  error  in  the  value  of  the  rotation  rate 
of  the  S-66  satellite  signals. 

In  this  study,  the  Tq  time  can  be  observed  experimentally  from  the 
satellite  passes.  The  accuracy  of  this  observation  was  limited  by  the 
width  of  the  recorder  marking  lines.  The  smallest  division  of  time  on 
the  marking  scale  was  5  seconds.  The  events  time  record  (see  Fig.  4) 
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FIG.  10  DIURNAL  VARIATION  OF  EQUIVALENT  SLAB  THICKNESS  AND  MEAN  SCALE  HEIGHT 
ASCENDING  AND  DESCENDING  PASSES 


can  be  resolved  to  within  1  second.  For  the  satellite  passes  with  no 

other  observable  T  time,  the  T  time  had  to  be  calculated.  When  cal- 
o  o 

culations  are  made  for  passes  with  observable  T  times,  the  calculated 

o 

values  of  T  are  almost  always  within  a  few  seconds  of  the  observed 
o 

values . 

The  total  errors  in  the  values  of  electron  content  derived  from  the 
rotation-rate  are  usually  less  than  10  percent.  However,  an  exception 
occurs  for  the  early  morning  hours,  when  the  rate  drops  to  less  than  1 
revolution  per  minute,  which  may  cause  errors  as  large  as  20  percent. 
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V  DEVELOPMENT  OF  THE  EQUATORIAL  ANOMALIES 


A.  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT 

The  latitudinal  variation  of  electron  content  is  determined  by  the 
total-rotations  method.  A  computer  program  has  been  prepared  for  Eq. 

(9)  that  is  used  to  compute  the  latitudinal  variation  of  electron  content 
for  the  selected  passes.  The  total  number  of  rotations  is  measured  from 
the  transverse  position  to  a  given  satellite  position,  and  the  corres¬ 
ponding  satellite  position  is  obtained  from  the  local-time  tabulations 
of  the  satellite  ephemeris.  From  this  calculation,  the  sub-ionospheric 
point  (the  projection  on  the  earth  of  the  intersection  of  the  ray  path 
between  satellite  and  ground  stations  with  the  centroid  of  the  iono¬ 
spheric  electron  distribution  profile)  traversed  about  10°  of  latitude 
around  a  ground  station.  Therefore,  the  latitudinal  variation  of  electron 
content  can  be  calculated  for  a  section  across  the  sky. 

To  increase  the  range  of  latitude  over  which  the  electron  content 
could  be  investigated,  two  remote  sites  were  set  up  at  Chiengmai  and 
Songkhla,  which  are  located  north  and  south  of  Bangkok  respectively 
(Fig.  11).  By  analyzing  data  collected  at  these  two  remote  stations 
as  well  as  at  Bangkok  from  November  1966  to  March  1967,  the  electron 
content  over  a  latitude  range  of  more  than  25°  could  be  obtained.  The 
selection  of  satellite  fading  records  that  would  yield  high  quality  for 
analysis  at  all  three  stations  was  governed  by  the  following  conditions: 

(1)  The  satellite  equatorial-crossing  longitude  for  ascending 
passes  should  be  between  98  and  107°E,  and  for  descending 
passes  should  be  between  93  and  102°E,  because  of  the 
inclination  of  the  satellite  orbit  and  the  locations  of 
three  stations. 

(2)  The  data-recording  period  must  be  long  enough  to  obtain 
data  overl aping  between  Songkhla  and  Bangkok,  and  between 
Chiengmai  and  Bangkok 

(3)  The  fading  records  should  be  free  of  scintillation. 
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FIG.  11  SHADED  ZONES  OF  SELECTED  SUB-IONOSPHERIC  PATH 
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The  electron  content  was  calculated  for  points  along  the  sub- 
ionospheric  path  (the  projection  on  the  earth  of  the  locus  of  the  point 
at  which  the  ray  path  between  satellite  and  ground  station  intersects 
the  centroid  of  the  ionospheric  electron-distribution  profile).  The 
sub-ionospher ic  paths  for  the  selected  passes  all  fall  within  the  shaded 
/.one  shown  in  Fig.  11.  The  /ones  for  ascending  and  descending  passes 
are  shaded  separately.  Records  for  which  the  electron  content  is  calcu¬ 
lated  were  selected  from  a  period  of  November  1966  through  March  1967. 

The  satellite  passes  occur  at  different  local  times  owing  to  the  pre¬ 
cession  of  the  orbit;  the  times  in  this  period  are  distributed  covering 
twenty-l'our  hours.  Therefore,  it  is  convenient  to  pool  the  selected 
passes  on  a  time-block  basis,  in  one-hour  groups.  Ephemeris  data  of  the 
selected  passes  for  each  hour  group  are  given  in  Appendix  C. 

The  results  of  the  latitudinal  electron  content  of  each  hour  group 
are  presented  in  Part  (a)  of  Figs.  12  through  32.  The  number  beside 
each  curve  is  a  key  to  the  data  on  which  the  satellite  pass  occurred; 
the  date  can  be  determined  by  reference  to  Appendix  D.  Notice  that 
the  data  for  the  local-time  blocks  0300-0400.  0500-0600,  and  1600-1700 
are  missing.  This  is  because  of  the  poor  quality  of  the  records  ar.d  the 
lack  of  observations  for  analysis.  The  computed  values  of  latitudinal 
electron  content  are  plotted  as  a  function  of  geographic  latitude  and 
magnetic  dip  angle.  The  principal  features  of  the  results  show  the 
diurnal  development  and  collapse  of  the  equatorial  anomalies  in  the 
electron  content.  The  latitudinal  variations  of  electron  content  from 
the  period  0000-0700  of  local  time  showed  a  rather  flat  variation  of 
latitudinal  distribution,  as  illustrated  in  Part  (a)  of  Figs.  12  through 
16.  These  figures  show  that  electron  content  at  the  southern  location 
is  higher  than  at  the  north.  This  effect  is  caused  by  the  season  when 
the  Faraday  fading  records  were  obtained,  since  the  records  were  observed 
for  the  period  of  winter  solstice.  The  sun  is  south  of  the  equator,  so 
the  expected  higher  values  of  electron  content  in  the  southern  hemisphere 
are  obtained.  At  the  period  0800-1000  hours  of  local  time,  the  latitudinal 
distribution  of  electron  content  is  characterized  by  a  peak  near  the 
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magnetic  dip  equator  as  illustrated  An  Part  (a)  of  Figs.  17  and  19.  The 
electron  content  variations  demonstrate  the  general  shape  of  the  equatorial- 
trough  anomaly  as  observed  between  1000  and  1600  local  time  [Part  (a)  of 
Figs.  20  through  25].  The  variations  show  a  trough  near  the  magnetic 
equator.  The  latitudinal  variation  of  electron  content  is  collapsed  to 
a  peak  and  shifted  north  of  the  magnetic  equator  near  the  sunset  hour 
1700-1800  local  time  as  illustrated  in  Part  (a)  of  Fig.  26,  and  the 
equatorial  trough  anomaly  is  developed  again  at  period  1800-0000  local 
time  as  illustrated  in  Part  (a)  of  Figs.  27  through  32.  The  exception 
of  one  peak  of  electron  content  near  the  magnetic  equator  as  shown  by 
Curve  1  in  Part  (a)  of  Fig.  29  is  not  well  understood.  This  diurnal 
development  and  collapse  of  the  equatorial  anomalies  has  been  reported 
by  other  workers.10  ,16 
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FIG.  13  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F-LAYER 
ORDINARY  CRITICAL  FREQUENCY,  0100-0200  LOCAL  TIME 
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FIG.  14  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F-LAYER 
ORDINARY  CRITICAL  FREQUENCY,  0200  -0300  LOCAL  TIME 
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FIG.  15  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F-LAYER 
ORDINARY  CRITICAL  FREQUENCY,  0400  -  0500  LOCAL  TIME 
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FIG.  17  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F-LAYER 
ORDINARY  CRITICAL  FREQUENCY,  0700-0800  LOCAL  TIME 
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FIG.  18  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F-LAYER 
ORDINARY  CRITICAL  FREQUENCY,  0800-0900  LOCAL  TIME 
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FIG.  19  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F-LAYER 
ORDINARY  CRITICAL  FREQUENCY,  0900-1000  LOCAL  TIME 
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FIG.  20  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F-LAYER 
ORDINARY  CRITICAL  FREQUENCY,  1000-1100  LOCAL  'I  I  ME 
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FREQUENCY  (toF2  )  -  MHz  ELECTRON  CONTENT - e*ecW»lO 
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FIG.  22  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F-LAYER 
ORDINARY  CRITICAL  FREQUENCY,  1200-1300  LOCAL  TIME 
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FIG.  23  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F-LAYER 
ORDINARY  CRITICAL  FREQUENCY,  1300-1400  LOCAL  TIME 
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FIG.  24  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F~LAYER 
ORDINARY  CRITICAL  FREQUENCY,  1400-1500  LOCAL  TIME 
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FIG.  26  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F-LAYER 
ORDINARY  CRITICAL  FREQUENCY,  1600-1700  LOCAL  TIME 
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FIG.  27  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F-LAYER 
ORDINARY  CRITICAL  FREQUENCY,  1800-1900  LOCAL  TIME 
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FIG.  29  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F-LAYER 
ORDINARY  CRITICAL  FREQUENCY,  2000-  2100  LOCAL  TIME 
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FIG.  30  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F~LAYER 
ORDINARY  CRITICAL  FREQUENCY,  2100  -  2200  LOCAL  TIME 
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FIG.  31  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F-LAYER 
ORDINARY  CRITICAL  FREQUENCY,  2200-2300  LOCAL  TIME 
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FIG.  32  LATITUDINAL  VARIATION  OF  ELECTRON  CONTENT  AND  F-LAYER 
ORDINARY  CRITICAL  FREQUENCY,  2300-0000  LOCAL  TIME 
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B.  LATITUDINAL  VARIATION  OF  F-LAYER  CRITICAL  FREQUENCY 


Data  collected  foA  calculating  electron  content  of  the  ionosphere  by 
the  Faraday  rotation  measurements  can  be  used  to  estimate  the  latitudinal 
variation  of  the  F-layer  ordinary-wave  critical  frequency,  foF2.  There 
have  been  many  attempts  to  show  that  the  model  ionosphere  used  for 
calculation  of  the  total  electron  content  is  not  too  critical — except  for 
the  refraction  correction.17-20  To  determine  the  latitudinal  variation 
of  the  F-xayer  ordinary  critical  frequency,  Chapman's  equation  was  chosen 
to  represent  the  electron  density  distribution  with  height  of  the  ionosphere. 
Although  several  very  questionable  assumptions  were  made  by  Chapman  in 
the  derivation  of  his  formula,  the  general  shape  of  the  electron  density 
distribution  with  altitude  as  given  by  Chapman's  equation  still  gives  a 

a  l  A  20  2 1 

reasonable  fit  to  the  actual  distribution  profile  as  we  know  it.  >  >  > 

A  simple  equation  has  been  proposed  by  Wright33  and  Potts14  for  the  iono¬ 
spheric  electron-density  distribution  profile.  This  equation  assumed 
that  the  shape  of  the  profile  was  similar  to  that  given  by  Chapman's 
equation  and  the  F-layer  ordinary  wave  frequency  (foF2),  and  can  be 
expressed  in  terms  of  total  electron  content  and  mean  scale  height  as 


f  oF2  = 


f 


Ndh 


,10 


|5 . 12  H  X  10 
Substituting  Eq.  (9)  into  Eq.  (17)  gives 


(17) 


f  oF2  = 


5.12  GH  X  10 


10 


(18) 


The  method  for  calculating  foF2  as  a  function  of  latitude  is  very 
similar  to  the  method  for  calculating  values  of  latitudinal  electron 
content.  A  computer  program  has  been  prepared  for  the  computation  of 
the  latitudinal  variations  of  foF2  values.  In  this  study,  the  value  of 
mean  scale  height  (H)  was  taken  from  the  data  of  Fig.  10,  which  includes 
the  diurnal  variation  of  H  at  Bangkok  derived  by  the  rotation-rate 
method  of  Faraday  analysis.  The  value  of  mean  scale  height  in  the  F 
region  of  the  ionospheric  profile  was  assumed  to  be  constant  in  latitude 
for  the  equatorial  region.  Allowing  this  present  assumption,  the 


52 


similarity  between  the  theoretical  and  experimental  values  of  mean  scale 
height  have  been  reported . 20 >33 ,24 

The  calculated  values  of  foF2  are  computed  from  the  3ame  selected 
Faraday  fading  records  that  were  used  for  the  analysis  of  latitudinal 
variation  of  total  electron  content.  The  results  are  shown  in  Part  (b) 
of  Figs.  12  through  32.  The  calculated  values  are  for  the  points  on  the 
sub-ionospheric  path  in  the  shaded  zones  of  Fig.  11.  The  latitudinal 
variations  of  foF2  showed  diurnal  development  and  collapse  of  the 
equatorial  anomaly.  This  anomaly  shows  that  the  magnetic  equator  is 
nearly  midway  in  the  latitude  range  depicted.  The  results  are  associ¬ 
ated  with  the  diurnal  development  and  collapse  of  the  equatorial  anomaly 
in  total  electron  content  value,  and  illustrate  the  similarity  between 
electron  content  and  foF2 .  The  results  are  represented  in  a  schematic 
illustration  in  Fig.  33,  which  shows  in  simplified  form  the  diurnal 
development  and  collapse  of  the  equatorial  anomalies  of  foF2  value  for 
the  winter  solstice.  During  the  period  0000-0700  local  time  the  lati¬ 
tude  variation  of  foF2  is  rather  linear,  with  the  southern-most  foF2 
values  significantly  higher,  as  illustrated  in  Part  (a)  of  Fig.  33.  The 
latitudinal  variation  of  foF2  is  characterized  by  a  peak  near  the  magnetic 
equator  during  the  period  0700-1000  local  time,  as  illustrated  in  Part  (b) 
of  Fig.  33.  The  well-defined  development  of  the  equatorial  anomaly 
occurred  during  the  period  1000-1600  local  time.  The  trough  near  the 
magnetic  equator  during  this  time  is  illustrated  in  Part  (c)  of  Fig.  33. 
The  equatorial  trough  anomaly  collapsed  during  the  next  interval.  Next, 
a  peak  formed  and  shifted  north  of  the  magnetic  equator  near  the  sunset 
hour  1700-1800  local  time,  as  illustrated  in  Part  (d).  The  equatorial 
trough  anomaly  is  developed  again  during  the  period  1800-0000  local  time, 
as  illustrated  in  Part  (e)  of  Fig.  33. 

To  check  the  validity  of  the  latitudinal  variation  of  foF2  as  deduced 
from  the  Faraday  rotation  data,  tne  results  were  compared  with  data  from 
vertical-incidence  sounders  taken  at  times  corresponding  to  satellite 
passages.  The  observed  foF2  values  at  Chiengmai  and  Songkhla  are 
obtained  from  Granger  sounders  which  were  operated  at  these  field  sites 
while  the  Faraday  measurements  were  made.  Due  to  the  lower  frequency 
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limit  of  these  sounders,  frequencies  below  4  MHz  are  not  available  for 
comparison.  A  comparison  of  the  calculated  and  observed  values  for  these 
three  stations  operating  ionosondes  is  presented  in  Appendix  D.  The 
calculated  values  are  for  the  points  on  the  sub-ionospheric  path  nearest 
the  ionosonde  locations  The  results  compare  well  with  the  observed 
values.  The  typical  difference  observed  is  less  than  10  percent. 
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VI  SUMMARY  AND  CONCLUSIONS 


This  report  presents  the  experimental  observations  of  the  equatorial 

ionosphere  made  from  Faraday  rotation  measurements.  Diurnal  variations 
of  electron  content  over  Bangkok  show  low  early-morning  values,  and 
maximum  values  in  the  afternoon.  Typically,  the  total  electron  content 

10  2  i 0  2 

varied  from  1  X  10  electrons/m  to  55  X  10  electrons/m  .  The  seasonal 
variations  show  the  maximum  values  of  electron  content  occurring  near 
the  equinoctial  periods  and  minimum  values  occurring  near  the  solstitial 
periods. 

The  results  of  the  latitudinal  variations  of  electron  content  and 
foF2  show  the  diurnal  development  and  collapse  of  the  equatorial  anoma¬ 
lies.  The  results  show  similarity  between  electron  content  and  foF2. 
During  the  period  0000-0700  local  time,  the  anomaly  in  latitudinal 
variation  is  rather  flat.  It  is  characterized  by  a  peak  near  the 
magnetic  equator  during  the  period  0700-1000  local  time.  A  well-defined 
development  of  the  equatorial  anomaly  occurs  during  the  period  1000- 
1600  local  time,  exhibiting  a  trough  near  the  magnetic  equator.  The 
equatorial  anomaly  collapses  to  a  peak  and  shifts  north  of  the  magnetic 
equator  near  the  sunset  hour  1700-1800  local  time;  the  equatorial-trough 
anomaly  develops  again  during  the  period  1800-0000  local  time. 

The  experimental  results  of  foF2  fro.n  Faraday  rotation  data  have 
been  compared  with  ionosonde  data,  usually  with  agreement  within  10 
percent.  The  latitudinal  variation  of  foF2  shows  the  diurnal  develop¬ 
ment  and  collapse  of  the  equatorial  anomaly.  The  scan  of  foF2  as  a 
function  of  latitude  can  be  used  to  adjust  the  foF2  observed  by  the  C-2 
ionosonde  at  Bangkok  for  the  prediction  of  maximum  useful  frequency 
(MUF)  in  other  parts  of  Thailand.  These  adjustments  are  needed  because, 
whf n  the  largest  gradients  occur  (between  the  development  and  collapse 
conditions  of  the  equatorial  anomalies),  the  variation  can  significantly 
affect  the  MUF  of  HF  radio  systems  in  Thailand.  It  has  been  shown  that 
observed  foF2  for  Bangkok,  when  available,  can  be  used  to  improve  local 
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frequency  prediction.*" ,2S  The  latitudinal  variations  of  foF2  value 
presented  provide  information  permitting  improved  predictions  of  MUF 
for  areas  ir  Thailand  where  ionosonde  data  are  not  available. 
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Appendix  A 

DERIVATION  OF  THE  GEOMETRICAL  COEFFICIENT  G-VALUES 
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Appendix  A 


DERIVATION  OF  THE  GEOMETRICAL  COEFFICIENT  G- VALUE 

The  method  of  calculation  of  the  G-value  is  given  in  this  appendix 
in  detail,  for  the  purpose  of  providing  a  computer  progr**..ii  for  direct 
application  to  the  analysis  of  the  Faraday  rotation  data  from  satellite 
radio  signals. 

The  G  value  is  calculated  at  an  ionospheric  centroid  height  and 
evaluating  the  magnetic  field  at  that  altitude.  Assume  no  refraction 
in  the  ionosphere  or  lower  atmosphere.  Consider  the  geometry  shown  in 
Figs.  A-l  and  A-2,  where  the  coordinates  of  the  subsatellite  point  A 


FIG.  A-l  GEOMAGNETIC-FIELD/RAY-PATH  GEOMETRY 
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FIG.  A- 2  SURFACE  ARC  DISTANCE 


(the  projection  of  satellite  on  the  earth),  the  subsatellite  point  E 

(the  projection  on  the  earth  of  the  intersection  of  the  ray  path  between 

satellite  and  ground  station,  with  the  centroid  of  the  ionospheric  electron 

density  profile),  and  the  ground  station,  R,  are  8g,  Xg;  8^,  X^;  and 

8  ,  X  ;  geographic  latitude  and  longitude,  respectively.  The  earth's 
r  r 

radius  is  denoted  by  R  and  the  ionospheric  centroid  height  is  h  . 

0  c 

Satellite  height  is  h  . 

s 

By  using  spherical  trigonometry,  the  angular  range  C  (measured  at 
the  center  of  the  earth)  is: 

cos  C  =  sin  8  sin  8  +  cos  8  cos  8  cos  (X  -  X  )  .  (A-l) 

s  r  s  r  s  r 

Let  the  slant  distance  from  satellite  to  ground  station  be  equal 

to  l: 

l  =  [R2  +  (R  +  h  )2  -  2K  (R  +  h  )  cos  C]?  .  (A-2) 

ees  ees 


62 


The  inclination  of  signal  ray,  I  ,  and  satellite  zenith  angle,  x>  at 

R 


the  ionospheric  centroid  height  can  be  determined  as 


IR  =  cot 


R  (R  +  h  )  sin  C 
e  e  s 


Ar  +h  )  [R  +(R  +h  )  -2R  (R  +h  )cos  Cl  -  [R  (R  +h 

VGC  g  0  g  60S  J  1  g  0 


)sln  C]2)1/2 


=  90°  -  I. 


(A- 3) 
(A-4) 


sec  x  =  cosec  I 


The  observing  satellite  zenith  angle,  Z,  is 


Z  =  90°  -  cos 


R  +  h 
e  c 


(A-6) 


Figure  A-2  shows  the  arcs  of  the  earth's  surface  from  the  geographical 
north  pole  to  the  subsatellite  point,  subionospheric,  point  and  ground 
stations.  (If  more  convenient,  we  may  join  them  to  the  geographical 
south  pole.)  Then  the  subionospheric  point  E  (9^,  is 

9J  =  sin  1  [sin  9  (cos  b  -  sin  b  cot  C)  +  sin  9  sin  b  cosec  Cl  (A-7) 
i  r  s 


\  =  \  +  N  (If  2  \  ) 

i  r  s  r 


(A-8) 


\ =  \  -  N  (If  \  S  X.) 

i  r  si 


where  f  t  1/2^ 

,  -1<  cln  b  f  1  5ln  a5  -  sln  9r  c°5  C1  1  > 

S  n  |^cos  I  cos  9^  sin  C  j  J  J 


(A-9) 


(A-10) 


b  -  cos 


R  (R  +  h  )  sin  C 

_ e  e _ s _ 

2  2  1  /2 
(R  +  h  )[R  +  (R  +  h  )  -  2R  (R  +  h  )  cos  C)  ' 

ece  es  ees 


(R  +  h  )  sin  C 

_ e  s _ 

2  2  1/2 
[R  +  (R  +  h  )  -  2R  (R  +  h  )  cos  C] 

e  e  s  ees 


(A-ll) 


The  declination  of  the  signal  ray,  D  ,  at  the  ionospheric  centroid 

R 

height  can  be  determined  by: 
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sin  8 

_ s _ 

cos  8  sin  (C  -  b) 

The  angle  between  signal  ray  path  and  earth's  magnetic  field,  8,  can 
be  determined  by  the  geometry  as  shown  in  Fig.  A-3.  Let  OX,  OY,  and  OZ 
be  three  perpendicular  lines  through  the  ionospheric  centroid  height  at 
Point  O,  where  I  and  are  magnetic  dip  angle  and  declination,  re¬ 
spectively.  Using  the  direction  cosines  theorem,  the  earth's  magnetic 
field  OB  is  in  the  direction  with  the  direction  cosines  (cos  a  .  cos  @  , 
cos  Yg),  and  the  signal  ray  OS  is  in  the  directions  (cos  cos 

cos  Y,,),  By  the  rule  for  the  addition  of  direction  cosines,  the  angle. 

8,  between  OB  and  OS  can  be  determined  by 

cos  8  =  cos  or  cos  a  +  cos  8  cos  3  +  cos  y  cos  y  .  (A-13) 

Ho  R  B  R  B 

By  the  spherical  trigonometry,  we  can  solve  for  the  cosine  value 
as  follows: 


cos 

aR  = 

cos 

XR 

cos 

°R 

cos 

0R  = 

cos 

XR 

sin 

dr 

cos 

33 

n 

sin 

XR 

(A-14) 

COS 

*B 

cos 

JB 

cos 

°B 

cos 

II 

CQ 

cn 

cos 

XB 

sin 

db 

cos 

YB 

sin 

XB 

• 

(A-15) 

Substituting  Eqs.  (A— 11)  and  (A-12)  into  (A-10), 

cos  8  =  cos  I  cos  I  cos  (D  -  D  )  +  sin  I_  sin  I  .  (A-16) 

R  B  R  B  R  B 

Note  that  I  ,  I  is  positive  when  the  direction  is  below  the 
R  B 

horizontal.  Since  FJg.  A-l  shows  that  I  (direction  of  signal  ray 

R 

lrom  satellite)  is  always  negative,  Eq.  (A-14)  becomes 

cos  8  =  cos  I  cos  I  cos  (D  -  D  )  -  sin  I  sin  I  .  (A-17) 

R  B  R  B  R  B 
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0* -4140-17* 


FIG.  A-3  RECTANGULAR  COORDINATE  OF  GEOMAGNETIC  FIELD/RAY  PATH 
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The  value  of  G  can  be  calculated  for  any  time  during  the  satellite 
pass  for  which  a  measurement  of  electron  content  Is  performed.  The 
expression  for  the  G  value  from  Eqs.  (A-5)  and  (A-17)  is  then 

G  =  B  cos  9  sec  ^  .  (A-18) 

The  value  of  B  in  Eq.  (A-18)  is  properly  specified  by  the  spherical 
harmonic  expansion  of  the  first  six  harmonics  to  which  about  0.3  percent 
of  the  local  geomagnetic  field  measurements  during  magnetically  quiet 
conditions.37  The  first  six  harmonics  of  spherical  harmonic  expansion 
are  described  in  detail  by  the  author  in  Ref.  28. 
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Appendix  B 

COMPUTER  PROGRAM  FLOWCHARTS 
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FIG.  B-2  SUBROUTINE  INPUT 
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Appendix  C 

SELECTED  SATELLITE  EPHEMERIS 
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Appendix  C 


Hour  Group 


0000-0100 

0100-0200 

0200-0300 


04  00-0500 


0600-0700 


0700-0800 


0800-0900 


The  letters  A 
respectively. 


15  January  1967 
22  January  1967 

20  January  1967 

7  January  1967 

8  January  1967 

10  January  1967 

4  October  1966 
20  December  1966 

5  December  1966 

8  December  1966 

9  December  1966 

11  December  1966 
13  December  1966 
19  December  1966 

5  March  1967 

29  November  1966 
7  December  1966 

2  March  1967 

3  March  1967 

6  March  1967 

7  March  1967 

28  November  1966 

30  November  1966 
1  December  1966 


and  D  indicate  ascending  and  descending  satellite  passes, 
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SELECTED  SATELLITE  EPHEMERIS  -  Continued 


Hour  Group 

Curve 

Number 

Date 

Revolut ion 

Number* 

Observing 

Time 

4 

22  February  1967 

11884D 

0850 

5 

25  February  1967 

11 92 5 A 

0820 

6 

26  February  1967 

11939A 

0850 

7 

28  February  1967 

11966A 

0800 

8 

1  March  1967 

11980A 

0830 

c> 

4  March  1967 

12021A 

0800 

0900-1000 

1 

15  February  1967 

11788A 

0900 

2 

16  February  1967 

11802A 

0930 

3 

20  February  1967 

11857A 

0930 

4 

23  February  1967 

11898A 

0910 

1000-1100 

1 

17  February  1967 

11816A 

1000 

1100-1200 

1 

4  February  1967 

11638A 

1105 

2 

5  February  1967 

11652A 

1130 

1200-1300 

1 

6  February  1967 

11666A 

1200 

1300-1400 

1 

15  January  1967 

11365A 

1410 

2 

20  January  1967 

11433A 

1300 

3 

21  January  1967 

11447A 

1330 

1400-1500 

1 

8  January  1967 

11269A 

1430 

1500-1600 

1 

10  January  1967 

11297A 

1625 

1700-1800 

1 

19  December  1966 

10997A 

1740 

1800-1900 

1 

13  December  1966 

10915A 

1830 

2 

14  December  1966 

10929A 

1850 

3 

20  December  1966 

1101 1 A 

1810 

| 

4 

6  March  1967 

12054D 

1825 

5 

9  March  1967 

12095D 

1810 

*The  letters  A  and  D  indicate  ascending  and  descending  satellite  passes, 
respectively. 
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SELECTED  SATELLITE  EPHEMERIS  -  Continued 


Hour  Group 


1900-2000 


Curve 

Number 


2000-2100 


2100-2200 


2200-2300 

2300-0000 


Date 

- , 

Revolut ion 

Number* 

1 - 

i  Observing 

Time 

4  December  1966 

10792A 

1930 

7  December  1966 

10833A 

1910 

8  December  1966 

10847A 

1930 

5  March  1967 

1 2  04 1 D 

1940 

29  November  1966 

1 0709A 

2020 

25  November  1966 

10668A 

2040 

1  December  1966 

10751A 

2000 

16  February  1967 

11808D 

2050 

2  March  1967 

12000D 

2010 

30  November  1966 

10737A 

2115 

13  February  1967 

11767D 

2110 

14  February  1967 

11781D 

2140 

17  February  1967 

11822D 

2110 

11  February  1967 

11740D 

2T00 

5  November  1967 

10395A 

2350 

*The  letters  A  and  D  indicate  ascending  and  descending  satellite  passes, 
respectively. 
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CALCULATED  AND  OBSERVED  F- LAYER  CRITICAL-FREQUENCY  (foF2)  VALUES 
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Obs  =  Observed  Values 

+  Observed  values  are  higher  than  calculated. 
_  observed  values  are  lower  than  calculated. 
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I  I  AUSTHA  ’T 

Para day  rotation  data  from  Satellite  S-66  were  accumulated  at  the  Electronics  Labora¬ 
tory  of  the  Military  Research  and  Development  Center,  Bangkok,  Thailand  from  November 
1964  through  June  1967,  a  period  of  intermediate  sunspot  number  (sunspot  minimum  con¬ 
dition  at  1964/1965  and  sunspot  maximum  condition  at  1968/1969).  These  data  provide  a 
meens  of  improving  ionospheric  frequency  predictions  for  regions  in  Thailand  remote 
from  Bangkok  (i.e.,  where  no  ionosonde  data  are  available). 

In  this  report,  two  methods  of  analysis — involving  rotation  rate  and  total  number  of 
rotations — are  applied  to  the  Bangkok  observations.  The  rotation-rate  method  is  used 
to  determine  the  electron  content  when  the  angle  between  the  ray  path  and  geomagnetic 
field  is  90*  (transverse  position).  The  total-rotations  method  is  used  to  determine 
latitudinal  variations  of  electron  content  of  the  equatorial  ionosphere  for  the  se¬ 
lected  passes  from  a  joint  analysis  of  observations  at  three  stations  in  Thailand-- 
Songkhla,  Bangkok,  and  Chiengmai. 

The  diurnal  variation  of  the  local  electron  content  shows  low  early-morning  values  and 
high  afternoon  values,  resulting  in  a  large  maximum-to-minimum  ratio.  A  secondary 
maximum  was  observed  during  the  nighttime.  The  diurnal  and  seasonal  variations  of 
electron  content  are  examined  together  with  the  variation  with  solar  activity  as 
measured  by  the  10.7-cm  solar  radio  flux. 

The  diurnal  development  and  collapse  of  the  equatorial  anomaly  in  electron  content  is 
clearly  seen  from  the  total-rotations  data.  The  Chapman  equation  is  used  for  conver¬ 
sion  of  these  total  electron  content  values  to  peak  electron  density  in  order  to 
compute  the  F-layer  critical  frequency  (foF2)  and  its  variation  as  a  continuous  func¬ 
tion  of  latitude.  The  computed  foF2  values  are  compared  with  those  obtained  from 
three  ionosondes  which  were  located  at  Songkhla,  Bangkok,  and  Chiengmai,  and  the 
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